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I .  INTRODUCTION 


Measured  temperature-dependent  quenching  coefficients  for  collislonal 

deactivation  of  HF(v^  *  I)  and  DF(v^  =  1)  molecules  by  HF(V2  ”  0)  or  DF(v2  * 

0)  molecules  are  available  because  these  systems  have  been  studied 

extensively There  is,  however,  a  scarcity  of  temperature-dependent 

quenching  coefficients  for  collislonal  deactivation  of  HF  species  and  HF 

isotopic  species  from  vibrational  states  above  the  first  vibrational  level. 

The  deactivation  from  the  upper  vibrational  levels  of  HF  and  DF  by  HF(v  *  0) 

or  DF(v  =  0)  at  room  temperature  has  been  studied  experimentally  by  direct 

laser  pumping  of  the  v  =  2  or  v  »  4  vibrational  levels  of  HF,2»^  low-pressure 

combustion  in  a  flow  cell,4’"*  reactive  flows  in  medium-pressure  flow  tubes,**-® 

o 

high-temperature  combustion  in  shock  tubes,  and  laser-induced  fluorescence  by 
sequential  photon  absorption.1®-1® 

In  a  recent  paper,  Wilkins  and  Kwok14  presented  a  kinetic  model  of 
infrared  laser-induced  fluorescence  experiments,  which  simulated  available 
experimental  data  for  vibrational  relaxations  of  HF(v^  *  1)  +  HF(v2  =  0)  and 
its  isotopic  analogs. 1^,16  ^at  rotational  nonequilibrium  model  was  based  on 
the  predicted  energy-transfer  mechanism  in  HF  systems  reported  in  a  trajectory 
study  by  Wilkins.12  The  model  explained  the  temperature  dependence  observed 
for  HF(v^  *  1)  vibrational  relaxation  by  HF(V2  *0).  In  fact,  the  entire 
temperature  dependence  of  the  reported  quenching  rate  coefficients  was 
explained  without  use  of  any  mechanisms  that  involve  dimerization. 

In  the  previous  trajectory  study12  on  HF(v^)  +  HF(v2  *  0)  collisions,  the 
following  conclusions  were  drawn: 


For  HF(v)  +  HF  collisions,  the  most  probable  paths  for  V  +  V 
exchange  involve  single  vibrational  quantum  transitions. 


2.  The  V  ■*■  R  processes  occur  by  means  of  an  intramolecular  energy- 
transfer  mechanism.  The  vibrationally  excited  HF  molecule  converts 
one  or  more  quanta  of  its  vibrational  energy  into  rotational 
energy,  thus  populating  its  high  rotational  states. 

3.  One  quenching  process  of  the  high  rotational  quantum  states  is  the 
slow  relaxation  by  R  +  (R',T*)  processes.  The  probability  of 
rotational  de-excitation  decreases  as  rotational  quantum  J 
increases.  At  high  values  of  J,  J  =  ±  1  rotational  transitions  are 
much  more  probable  than  multiple  quantum  rotational  transitions. 

4.  HF  quasidimers  formed  at  room  temperature  or  above  are  not 
necessary  to  account  for  the  fast  V  +  R  self-relaxation  rates  with 
inverse  temperature  dependences  measured  by  various  experimental 
techniques. 

Additional  trajectory  studies^*^  verified  that  the  conclusions  reached 
for  the  mechanisms  in  HF(vj)  +  HF(v2)  collisions  are  equally  applicable  to 
DF(v^)  +  DF(v2)  and  HF(v^)  +  DF(v2)  systems.  The  exception  to  these 
mechanisms  involves  the  lack  of  V  -*■  V  energy  transfer  in  DF(v^)  +  HF(v2  =  0) 
systems.  However,  this  result  has  been  confirmed  in  laser-excited  induced 
fluorescence  experiments.  The  large  endothermicities  involved  in  DF(vi)  + 
HF(v2  =  0)  collisions  probably  prevent  the  occurrence  of  the  V  ♦  V  process. 

The  collision  dynamics  indicate  that  in  V  +  V  transfer  of  DF(v^)  +  HF(v2>  and 
HF(v^)  +  DF(v2)  collisions,  vibrationally  excited  HF(v^)  molecules  transfer 
vibrational  energy  to  the  DF  molecules. 


The  previous  three  trajectory  studies  reached  the  following  conclusions 
concerning  V  -*■  V  transfer: 

1.  In  both  HF(v^)  +  HF(v2  =  0)  and  DF(v^)  +  DF(V2  =  0)  systems,  the 
rate  coefficients  for  the  V  V  processes  decreased  as  v  increased 
because  of  the  increased  endothermlcities  with  increasing  v. 

2.  In  HF(v^)  +  DF(v2  *  0)  collisions,  the  rate  coefficients  for  V  ♦  V 
processes  increased  slightly  as  v  increased  because  of  the  decrease 
in  exothemicities  with  increasing  v. 

3.  In  DF(v^)  +  HF(V2  “  0)  collisions,  rate  coefficients  for  V  ♦  V 
processes  can  be  neglected,  since  they  are  calculated  from  small 
cross  sections  that  result  from  their  large  endothermlcities. 

4.  The  rate  coefficients  for  the  V  +  V  processes  had  a  temperature 
dependence. 

It  is  difficult  to  assess  separate  actual  values  for  V  V  or  V  R 
processes  from  the  experimental  measurements,  because  measured  coefficients 
usually  represent  the  sum  of  the  contributions  from  the  V  +  V  and  the  V  R 
energy- transfer  processes  for  many  transfer  channels. 

In  addition,  the  calculated*^"*^  temperature-dependent  rate  coefficients 
cannot  be  compared  directly  with  the  experimental  temperature-dependent  rate 
coefficients  for  the  vibrational  deactivation  of  HF(v^)  or  DF(v^)  by  either 
HF(v2  =  0)  or  DF(v2  *  0) .  The  latter  should  be  called  "empirical  quenching 
coefficients,"  since  they  represent  net  removal  rates  of  HF(vj)  or  DF(v^) 
through  a  number  of  energy-transfer  channels.  Previous  experimental  analyses 
of  deactivation  of  these  species  from  upper  vibrational  states  did  not 
properly  consider  the  role  of  V  -*•  R  and  R  ♦  V  energy-transfer  processes. 
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The  latter  mechanisms  are  believed  to  be  important  in  the  vibrational 
deactivation  of  hydrogen  halides,  particularly  HF  and  its  isotopic  species, 
from  upper  vibrational  levels.  Most  current  experimental  methods  are  not 
capable  of  directly  measuring  the  fine  structure  rate  coefficients  for  the 
energy-transfer  processes  as  predicted  by  the  trajectory  studies. 

The  rate  coefficients  for  V  ♦  R,  R  *  V,  R  (T',R')  and  (R',T")  *  R  and 
V  +  V'  energy-transfer  processes  obtained  from  the  trajectory  studies  were 
incorporated  into  a  nonequilibrium  kinetic  computer  program  that  models  the 
laser-induced  fluorescence  technique.  The  purpose  of  this  study  is  to 
determine  if  this  rotational  nonequilibrium  model  can  reproduce  the 
experimental  temperature-dependent  quenching  rate  coefficients  for  vibrational 
deactivation  of  HF(v^  =  2  to  6)  by  HF(v2  =  0),  HF(vj  =  2  to  6)  by  DF(v2  =  0), 
DF(v^  =  2  to  6)  by  DF(v2  =  0),  and  DF(v^  =  2  to  6)  by  HF(v2  =  0).  The 
quenching  rate  coefficients  for  the  V  ♦  R  processes  are  obtained  by 
subtracting  the  rate  coefficients  for  the  approximate  V  V  processes  from  the 
total  decay  rate  (i.e.,  kq.  =  kw  +  .  The  deactivation  rates  for  the  V  +  R 

processes  are  then  used  to  determine  how  V  +  R  rates  scale  with  vibrational 
level  v. 
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II.  MODEL 


In  a  previous  paper,  Wilkins  and  Kwok^  described  the  computer  code 
developed  to  calculate  time-dependent  HF(v^,  J^)  concentrations  in  defined 
vibrational  (vp  and  rotational  (Jj)  states  for  specific  kinetic  experiments 
relevant  to  energy  transfer  in  HF.  The  description  of  the  model  presented  in 
that  paper  will  not  be  repeated  here.  The  model  is  extended  to  higher  v^ 
levels  by  incorporation  of  the  additional  fine  structure  rate  coefficients  for 
V  -*•  R,  R  V,  R  +  (R',T')  and  (R',T')  ■>  R  and  V  *  V’  obtained  from  the  tra¬ 
jectory  studies. 15-17  t^e  appencjices  are  tabulated  rate  coefficients  em¬ 
ployed  in  the  model  calculations.  The  rate  coefficients  for  rotational  relax- 
ation  are  assumed  identical  in  each  of  the  v^  manifolds.  The  NEST  code1  is 
used  to  calculate  the  rotational  nonequilibrium  effects.  Time-dependent 
HF(vj,  Jj)  densities  must  be  generated  by  this  model  for  comparison  with 
actual  experiments.  The  computer  number  densities  are  summed  over  J,  and  the 
sum  is  analyzed  in  the  same  manner  as  in  actual  experiments  to  deduce  the  em¬ 
pirical  quenching  coefficients  reported  by  experimenters  (generally  for  room 
temperature) . 

The  computer  model  employed  to  obtain  total  decay  rates  for  relaxation 
from  the  upper  vibrational  levels  of  HF(vp  simulates  a  novel  variation  of  the 
method  of  laser-induced  fluorescence. ^  The  experimental  technique  involves 
exciting  the  first,  the  first  and  second,  the  first,  second,  and  third,  etc. 
vibrational  levels  of  HF  with  the  use  of  multiband  laser  optical  pumping.  The 
vibrational  level  of  interest  is  then  populated  by  collisions  of  the  photo- 
lytically  excited  HF  species.  The  mathematical  model  presented  by  Osgood, 
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Sackett,  and  Javan*®  is  used  here.  They  demonstrated  that  the  temporal 
evolution  of  ny(t)  can  be  written  as 

ny(t)  -  A(  1  -  e‘t/TR)  e^l*  (1) 

where  A  is  a  constant,  xR  is  the  risetime,  and  is  the  decay  rate  from  the 
v^  *  1  level.  The  temporal  behavior  of  nv  is  described  by  a  factor  with  a 
rising  exponential  whose  time  constant  is  essentially  the  total  decay  rate  for 
de-excitation  from  level  v  times  a  factor  with  an  exponential  decay  rate  v 
times  the  V  ■*  R  decay  rate  from  the  v  <■  1  level.  xD  is  defined  by 

^  "  kw  +  kvR  -  S  (2> 


where 


k 


T 


k 

w 


+  k 


vR 


(3) 


and  kj.  is  the  total  decay  rate  from  level  v,  kvy  is  the  contribution  to  kj. 
from  V  +  V  energy-transfer  processes,  and  k^  is  the  contribution  to  kq-  from  V+ 
R  processes.  To  obtain  xR,  the  time  at  which  the  maximum  fluorescence 
intensity  occurs  is  measured.  The  equivalent  is  obtained  theoretically  by 
solving  for  the  point  in  time  at  which  i^Ct)  is  a  maximum.  Setting  dn^tj/dt 
equal  to  zero,  and  solving  for  tmax  we  obtain 

/  xH(v>\ 

t  “  x_ln  (l  + - )  (4) 

max  R  \ 
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where  the  decay  time  t.(v)  -  1/vy,  .  The  risetime  t_  is  obtained  by  solving 
the  transcendental  equations  numerically.  This  technique  is  applicable  as 
long  as  tr  <  <  Td. 

The  initial  conditions  for  each  computation  depend  on  proper  adjustment 
of  the  laser-excited  population  of  the  levels  n°.  Osgood,  Sackett,  and 
Javan*®  demonstrated  that  the  time -dependent  amplitudes  of  ny  at  vibrational 
equilibrium  are  defined  by 

2b)  X  v 

o  e  e  v-"  .  . 

*  °1  6Xp  ~ kT~  ^  (j  " 

[  j-i 

where  o^x^  is  the  anharmonicity  constant  for  the  diatomic  molecule.  This 

choice  of  n°  removes  all  dependence  on  the  previously  measured  kvv  +  k^  and 

decouples  the  behavior  of  the  vth  level  from  that  of  all  higher  vibrational 

levels.  At  a  given  temperature,  pressures  of  HF(V2  »  0)  ■  H  [HF(v2  *  0,J2>] 

J 

up  to  0.5  Torr  are  used.  In  the  nomenclature  of  Osgood,  Sackett,  and 
Javan,  n°/nQ  was  taken  to  be  0.01  at  temperatures  from  300  to  700  K,  0.1  from 
700  to  1400  K,  and  0.2  from  1600  to  2400  K. 

As  in  our  previous  simulation,*^  extreme  care  has  been  taken  to  duplicate 
the  reported  initial  conditions  of  the  experiments  [Eq.  (5)]  as  well  as  the 
interpretive  approaches  of  the  experimenters  [Eqs.  ( 1 ) — ( 4) ]  as  we  relate  model 
number  density  results  to  their  works.  Then  we  extend  the  temperature  range 
beyond  that  of  their  experiments. 


III.  RESULTS  AND  DISCUSSION 


As  with  the  actual  experiments,  the  rise  rate  xR  is  obtained  from  Eq. 

(4).  The  calculation  is  performed  at  several  HF  densities  to  verify  that  the 
total  decay  rate  (kvv  +  ky^)  is  the  result  of  HF  quenching.  The  slope  of  the 
inverse  characteristic  decay  times  versus  HF(v  *  0)  densities  yields  the  same 
computed  value  of  kyv  +  kyR.  This  process  has  been  used  by  experimenters  to 
ensure  that  the  deduced  observed  quenching  coefficient  k^  +  kyR  is  not 
affected  by  competing  secondary  processes. 

Table  I  presents  a  summary  of  the  number  of  species  HF(v^,  J^)  and  rate 
coefficients  for  the  different  energy-transfer  processes  involved  in  the  use 
of  this  method  to  compute  the  total  decay  rate  from  the  v  =  3  level  of  HF.  A 
plot  of  a  typical  output  from  the  computer  for  tiy(t)  with  v  «  1,  2,  and  3  as  a 
function  of  time  is  given  in  Fig.  1.  From  Fig.  1,  t^  *  5.0  x  10-6  sec, 

A  — i  ->4  —1 

xd(v  *  1)  »  Yj  »  2.5  x  10  sec  A,  x^(v  =  2)  =  5.0  x  10  sec  ,  and  x^v  «*  3)  ** 
7.5  x  10^  sec-*.  It  was  found  that  x^(v  *  2)  and  T^(v  **  3)  were 
approximately  2x^(v  =  1)  and  3x^(v  *  1),  respectively.  Application  of  Eq. 
(4)  yields  xR  *  2.75  x  10  ^  sec  and  kyV  +  kyR  ■  17.5  x  10*^  cm^/mole  sec. 

From  the  trajectory  study, kyy  for  v  ■  3  at  300  K  has  a  value  of  5.2  x  10*^ 
cnr/mole  sec.  If  we  subtract  kyy  from  kj.,  we  obtain  kyR  *  12.3  x  10 
cm^/mole  sec  at  room  temperature.  The  inverse  of  the  quenching  coefficients 
for  V  ♦  R  processes  for  HF(v)  +  HF  and  its  isotopic  analogs  is  illustrated  on 
traditional  Landau-Teller  plots  in  Figs.  2  through  5,  where  Px(v)  is  plotted 
as  a  function  of  T-*^.  The  model  described  by  Wilkins  and  Kwok*^  replicated 
the  empirical  quenching  coefficients  for  HF(v  ■  1)  +  HF  relaxation.  The  model 
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Total  decay  rate  from  the  v  *  3  level  of  HF.  Number  density 
n^t)  versus  t(ysec).  The  first  and  second  vibrational  levels  of 
HF  are  excited  by  multiband  laser  optical  pumping.  The  v  *  3 
level  is  populated  by  collisions  of  the  photolytically  excited  HF 
species. 
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Fig.  3.  Rate  coefficients  for  vibrational  relaxation  of  DF(vj  «■  1 

through  6)  by  DF(vo  ■  0).  PT(iiaec-atm)  versus  T  and  T-*' ^  with  T 
in  K. 
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4.  Rate  coefficients  for  vibrational  relaxation  of  HF(v,  -  1 
through  6)  by  DF(v2  *  0).  Pt (psec-atm)  versus  T  and  T~*' ^ 
with  T  in  K. 
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Fig.  5.  Rate  coefficients  for  vibrational  relaxation  of  DF(v,  -  1 
through  6)  by  HF(v2  ”  0).  Px(Msec-atm)  versus  T  and 
with  T  in  K. 


used  in  this  paper  predicts  a  temperature  dependence  for  relaxation  from  the 
upper  vibrational  levels  similar  to  that  obtained  for  relaxation  from  the  v  ■ 

1  level.  There  are  insufficient  experimental  data  to  substantiate  this 
prediction  of  temperature  dependence  for  vibrational  relaxation  from  the  upper 
vibrational  levels.  A  comparison  of  the  ratios  ltvR(v)/^VR(v  ”  1)  f°r  the  V  + 

R  processes  indicates  that  the  rate  coefficients  for  deactivation  from  the 
upper  v  levels  scale  as  v11  with  n  varying  from  2.3  to  1.6  as  v  varies  from  2 
to  6.  The  actual  values  of  n  obtained  from  this  model  for  the  four  isotopic 
systems  of  HF(vj)  +  HF(v2  ■  0)  are  given  in  Table  II.  The  index  n  is  not  a 
constant  for  all  v.  The  model  predicts  only  a  slight  increase  in  rate 
coefficients  for  V  -*•  R  relaxation  as  v  increases  beyond  v  **  5.  This  model 
indicates  that  the  exponent  n  is  not  a  constant,  but  decreases  as  v 
increases.  The  experimental  values  of  n  are  obtained  from  total  decay 
quenching  coefficients  which  include  contributions  from  both  V  +  V  and  V  R 
processes.  The  vn  scaling  employed  in  our  model  refers  only  to  V  +  R 
processes.  The  V  +  V  contribution  can  decrease  or  increase  with  increasing  v, 
depending  mainly  on  the  behavior  of  the  energy  defects  with  increasing  v.  The 
vn  scaling  for  V  R  processes  is  expected  to  exist  even  in  systems  where 
initially,  at  low  v,  the  contribution  from  V  -*•  R  processes  is  extremely  small 
compared  with  the  contributions  from  V  -*•  V  processes.  Table  III  provides  a 
comparison  of  total  deactivation  rates  at  room  temperature  for  the  HF(vj)  + 
HF(v2  *  0)  system.  The  results  of  Kwok  and  Wilkins^  for  v  »  2  to  3  are  in 
good  agreement  with  those  of  Airey  and  Smith^*  and  Poole  and  Smith, while  the 
results  of  Osgood,  Sackett,  and  Javan*®  are  a  factor  of  2  larger.  The  v  -  4 
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TABLE  II.  V-dependence  of  V  +  R  rate  coefficients,  values  of  superscript  "n 


f 


,  V  +  R 
k 


V  +  R 
1+0 


Systems 


HF(v)  +  HF(0)  DF(v)  +  HF(0)  HF(v)  +  DF(0)  DF(v)  +  DF(0) 


TABLE  III.  Comparison  of  quenching  coefficients  at  T  *  300  K 
for  the  HF(v^)  +  HF  (v£  *  0)  system3 


energy-transfer  processes. 


level  can  be  considered  pivotal  for  determining  a  scaling  trend  in  v  ,  since  V+ 
R  processes  are  considered  dominant  over  V  -*•  V  processes  (83%  V  R  processes 


according  to  Table  IV).  However,  the  experimental  results^ vary  over 

nearly  a  factor  of  3.  With  higher  v,  any  experimental  scaling  trend  should  be 

still  more  definitive,  but  results  for  a  given  v  differ  by  a  factor  of 
2>4, 5, 7, 22 

The  model  used  here  gives  HF(v^)  +  HF(v2  =  0)  results  which  are  in  good 

agreement  with  the  flow  cell  experimental  measurements  of  Poole  and  Smith, ^ 

even  to  the  trend  of  exponent  n  with  increasing  v.  The  model  agrees 

qualitatively  with  Osgood^®  at  v  =  2,  3  but  not  at  all  at  v  =  4.  More  recent 
3  5  7  22 

results  *  *  »  indicate  that  a  coupling  problem  may  exist  on  Osgood's  v  =  4 

measurements.  The  results  of  Kwok  and  Wilkins^  agree  with  the  model,  but 

23  4  27 

those  of  Kwok  and  Cohen  do  not.  Both  Airey  and  Kwok  imply  quantitatively 

OO 

larger  quenching  coefficients,  and  the  results  of  Kwok  provide  a  nonvarying 
n  with  increasing  v. 

The  model  and  virtually  all  experiments  agree  that  the  exponent  n  is 

approximately  2  and  not  unity.  This  value  indicates  that  the  quenching  is  a 

complex  manifestation  of  several  mechanisms  and,  therefore,  cannot  be 

understood  as  single  quantum  transfer  between  simple  harmonic  oscillators. 

The  model  framework  demonstrates  that  this  complexity  in  quenching  coefficient 

behavior  may  be  the  integrated  effect  of  many  transfer  channels  and 

multiquantum  sets  of  V  *  R,  V  ♦  V,  and  R  R,  T  mechanisms.  Recently, 

3 

Douglas  and  Moore  measured  vibrational  relaxation  of  HF(v  =  3  and  4)  by 
pumping  HF  directly  from  v  =  0  to  v  =  4.  Table  III  demonstrates  that  their 
value  for  the  total  decay  rate  from  v  *  4  is  almost  two  and  one-half  times  as 
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TABLE  IV.  Percent  Contribution  of  v  +  R  energy  transfer  to  the  sum  of 
the  (v  ♦  v  +  v  R)  relaxation  in  HF(v^)  +  HF(v2  *  0) 
collision 


V 

300  K 

500  K 

700  K 

1000 

2 

37.0 

42.6 

45.2 

48.5 

3 

70.3 

73.9 

73.5 

76.2 

4 

83.4 

86.0 

85.3 

87.9 

5 

92.4 

93.9 

93.6 

94.3 

6 


96.1 


96.7 


96.5 


96.8 


large  as  its  value  from  v  »  3.  Douglas  and  Moore  found  no  evidence  for 
multiquantum  transitions  in  their  study  of  vibrational  relaxation  of  HF(v  *  3, 
4).  They  deduced  that  75  to  95 %  of  the  HF(v  =  4)  must  appear  as  HF(v  ■  3) 
by  observing  unresolved  low  J  states  and  assuming  rotational  equilibrium.  It 
is  not  clear  if  one  can  determine  that  the  quenching  process  is  a  single 
vibrational  quantum  process  by  observing  only  the  low  J-state  densities  of 
HF  in  only  the  v  *  3  and  4  levels. 

The  model  results  can  be  compared  with  available  theoretical 
models.  1^~21  In  tj,e  shin  and  Kim^  model,  it  is  assumed  that  only  V  +  V 
transfer  is  important  and  that  V  -*■  V  transfer  occurs  in  single-quantum 
transitions.  At  low  temperatures,  the  model  assumes  a  nonrigid  dimer 
formation  and  at  higher  temperatures,  complete  rotational  motion  of  the 
colliding  molecules.  The  energy  defect  is  removed  by  this  rapid  rotation. 

The  total  decay  rate  is  the  sum  of  the  two  mechanisms.  The  model  used  in  this 
paper  disagrees  with  the  Shin  and  Kim  predictions  and  predicts  that  V  R 
processes  make  considerable  contributions  to  the  total  decay  rates. 

Billing  and  Poulsen^0  used  a  semiclassical  approach  to  calculate  rate 
coefficients  for  vibrational  relaxation  of  HF(v)  +  HF(O)  with  v  =  1  through 
7.  They  treated  the  rotation  and  translation  classically  and  vibrational 
energy  transfer  quantum  mechanically.  They  found  that  orbiting  collisions 
were  important.  Three-dimensional  trajectory  studies  performed  by  Wilkins  do 
not  support  their  conclusion  concerning  the  role  of  orbiting  collisions  at 
room  temperature  and  above.  Even  though  their  model  is  different  from  our 
model,  which  predicts  formation  of  high  rotational  states  in  the  v  manifolds, 
their  model  did  predict  that  V  *  R  processes  make  a  dominant  contribution  to 
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/ 


T 


the  total  decay  rates.  Their  absolute  values  of  total  decay  rates  were  about 

55Z  of  the  experimental  values.  However,  It  Is  Important  for  the  production 

of  high  rotational  states  to  be  Investigated.  Such  states  might  Indicate  the 

validity  of  multiquantum  vibration-rotation  energy  conversion.  The  model  used 

20 

by  Billing  and  Poulsen  predicts  single-quantum  V  -*■  R  energy-transfer. 

21 

Clendening  et  al.  used  surprisal  theory  to  calculate  total  decay  rates 

for  HF(v)  +  HF(0).  This  method  allows  them  to  Include  contributions  from  rate 

coefficients  for  several  channels.  The  total  deactivation  rates  agree  well 

with  the  experimental  data  in  Table  III.  Their  total  decay  rates,  however, 

scale  linearly  with  v  in  the  same  way  as  an  harmonic  oscillator.  The  model 

presented  here  for  V  -*■  R  processes  scales  more  rapidly  with  v.  The  Billing 
20 

and  Poulsen  model  also  indicates  more  rapid  scaling  with  v  than  that 
indicated  by  Clendening  et  al. 

Table  IV  lists  the  percentage  of  the  total  decay  rates  for  HF(v)  +  HF(0) 
that  corresponds  to  V  R  energy  transfer.  At  300  K,  this  percentage 
increases  from  37%  at  v  *  2  to  96%  at  v  *  6.  Since  our  V  *  V  rates  decrease 
with  increasing  temperature,  the  percentage  contribution  that  results  from 
V  -*•  R  processes  also  increases  with  increasing  temperature  at  a  specific  v 
level.  At  1000  K,  the  percentage  contribution  to  the  total  decay  rate  that 
results  from  V  +  R  processes  increases  from  49%  at  v  «  2  to  97%  at  v  *  6. 

These  results  have  not  been  verified  experimentally.  Table  IV  indicates  that 
as  v  increases,  the  main  contribution  to  the  total  decay  rate  comes  from  V  R 
processes  and  not  from  the  V  ♦  V  processes.  This  conclusion  was  reached 
experimentally  by  Rwok  and  Wilkins^  and  Poole  and  Smith^  while  analyzing  their 
flow-tube  experiments,  and  theoretically  by  Billing  and  Poulsen.^®  Table  IV 
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indicates  that  at  T  *  300  K  and  v  ■  2,  the  contribution  from  V  V  processes 
is  approximately  63%  of  the  sura  +  kvR.  With  increasing  v,  however,  the 
V  +  R  processes  dominate.  A  previous  trajectory  study  by  Wilkins^  predicted 
that  the  V  ♦  V  rate  coefficients  for  HF(v)  +  HF(v  *  0)  collisions  decrease  as 
v  increases.  The  net  result  is  that  both  theoretical  and  experimental  deter¬ 
minations  for  HF(v)  +  HF(0)  collisions  agree  that  the  V  +  V  contributions  to 
the  total  decay  rate  coefficients  decrease  with  increasing  v. 

Table  V  provides  a  comparison  of  the  total  decay  rates  obtained  from  this 

model  with  available  experimental  data  for  the  HF  +  DF,  DF  +  DF,  and  DF  +  HF 

22 

systems.  The  values  measured  by  Kwok  for  HF(v)  +  DF  at  v  *  5  and  6  are  a 
factor  of  2  larger  than  the  model  predictions.  For  the  DF(v)  +  DF  system  at  v 
=  2,  the  comparison  between  this  model  and  experiment  is  good.  At  low  v  for 
the  DF(v)  +  HF  system,  the  model  predictions  are  in  good  agreement  with 
experimental  results  of  Bott  and  Kwok  and  Wilkins.  At  high  v  the  model 
predictions  are  almost  a  factor  of  2  larger  than  the  measurements  of  Poole  and 
Smith. To  the  best  of  our  knowledge,  Table  V  gives  all  the  experimental  data 
available  on  deactivation  rates  from  the  upper  v  levels  for  the  three  systems 
HF  +  DF,  DF  +  DF,  and  DF  +  HF.  The  model  has  been  used  to  calculate 
temperature-dependent  rate  coefficients  for  the  V  •*•  R  processes  for  these 
three  systems;  the  results  are  given  in  Figs.  3  through  5.  These  results 
should  be  very  useful  for  the  modeling  of  either  HF  or  DF  chemical  lasers. 
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TABLE  V.  Comparison  of  HF-DF,  DF-DF  and  DF-HF  total 
relaxation  rates  at  T  “  300  K 


IV.  CONCLUSIONS 


The  room  temperature  quenching  rate  coefficients  for  relaxation  from  the 
upper  vibrational  levels  of  HF  and  DF  have  been  successfully  duplicated  with 
good  agreement  with  available  experimental  data.  This  agreement  was  obtained 
with  the  use  of  a  rotational  nonequilibrium  model  and  rate.coef f icients  com¬ 
puted  by  Wilkins.  ^-17  ^ey  processes  are  V  /*•  R  and  R  *  V  mechanisms, 

which  give  the  problem  a  multiple-channel  nature.  The  model  provided 
temperature-dependent  quenching  rate  coefficients  for  relaxation  of  HF(v)  and 
DF(v)  by  HF(v  *  0)  and  DF(v  =  0).  The  temperature  dependence  is  predicted  to 
be  similar  to  that  measured  for  HF(v  =  1)  +  HF  and  its  isotopic  analogs.  The 
results  of  this  study  indicate  that  HF(v)  +  HF(v  =  0),  DF(v)  +  DF(v  =  0),  ^ 


DF(v)  +  HF(v  *  0),  HF(v)  +  DF(v  *  0)  all  scale  as  vn  with  n  varying  between 
2.3  and  1.6  as  v  increases  from  2  to  6. >  This  scaling  has  not  been  verified 
experimentally  for  temperatures  higher  than  room  temperature.  The  scaling  is 
applicable  only  to  the  V  ♦  R  processes.  For  example,  in  HF(v)  +  DF(v  •  0), 
the  V  +  V  process  was  predicted  to  increase  slightly  with  increasing  v,  but 
the  V  R  processes  were  found  to  scale  similarly  to  those  for  HF(v)  +  HF(v  = 
0).  With  the  appropriate  fine  structure  rate  coefficients,  this  model  should 
be  equally  applicable  to  vibrational  relaxations  from  the  upper  vibrational 
levels  of  other  hydrogen-halide  molecules.  This  theoretical  study  is  the 
first  in  which  the  temperature  dependence  of  the  V  +  R  rate  coefficients  for 
HF(vj)  +  HF(v2  *  6)  and  its  isotopic  analogs  has  been  predicted  over  the 
entire  temperature  range  from  300  to  2400  K. 


For  a  given  molecule  such  as  HF(v),  the  insensitivity  of  quenching 
coefficient  to  value  of  the  exponent  n  for  a  number  of  differing  chaperones, 
the  absolute  value  of  n  itself  (near  2),  and  the  insensitivity  of  the  relative 
temperature  dependence  of  the  quenching  coefficient  regardless  of  v  level  all 
indicate  that  the  previously  observed  quenching  coefficient  is  not  sensitive 
to  details  of  process  cross  sections.  There  is  a  commonality  in  the  details 
of  the  internal  energy  structure  of  the  molecule.  This  Implies  that  much  of 
this  type  of  energy  exchange  is  intramolecular  and  that  a  large  number  of 
channels  may  be  involved.  Certainly,  questions  about  the  utility  of  the 
quenching  coefficient  in  revealing  the  actual  energy  transfer  processes  can  be 
raised. 


32 


REFERENCES 


(a)  J.  F.  Bott,  J.  Chem.  Phys.  J57,  96  (1972);  61,  3414  (1974). 

(b)  J.  J.  Hlnchen,  J.  Chem.  Phys.  _59,  2224  (1973);  _59,  233  (1973). 

(c)  J.  F.  Bott  and  N.  Cohen,  J.  Chem.  Phys.  55,  3698  (1971). 

(d)  W.  C.  Solomon,  J.  A.  Blauer,  F.  C.  Jaye,  and  J.  G.  Hnat,  Int.  J. 

Chem.  Klnet,  3_,  215  (1971). 

(e)  L.  S.  Blair,  W.  D.  Breshears,  and  G.  L.  Schott,  J.  Chem.  Phys.  59, 
1582  (1973). 

(f)  R.  A.  Lucht  and  T.  A.  Cool,  J.  Chem.  Phys.  ^3,  3962  (1975);  j>0,  1026 
(1974);  60,  2554  (1974). 

(g)  R.  R.  Stephens  and  T.  A.  Cool,  J.  Chem.  Phys.  56,  5863  (1972). 

(h)  J.  L.  Ahl  and  T.  A.  Cool,  J.  Chem.  Phys.  58,  5540  (1973). 

(1)  J.  R.  Airey  and  S.  F.  Fried,  Chem.  Phys.  Lett.  _8,  23  (1971). 

(j)  S.  S.  Fried,  J.  Wilson,  and  R.  L.  Taylor,  IEEE  J.  Quantun  Electron. 
QE-9,  59  (1973). 

(k)  R.  M.  Osgood,  P.  B.  Sackett,  and  A.  Javan,  J.  Chem.  Phys.  60,  1464 
(1974). 

(l)  M.  A.  Kwok  and  R.  L.  Wilkins,  J.  Chem.  Phys.  63  ,  2453  (1975). 

(m)  J.  K.  Hancock  and  W.  H.  Grier,  J.  Chem.  Phys.  57,  4515  (1972);  56, 
2474  (1972). 

M.  J.  Blna  and  C.  R.  Jones,  Appl.  Phys.  Lett.  22 ,  44  (1973). 

J.  Douglas  and  C.  B.  Moore,  Chem.  Phys.  Lett.  57,  455  (1978). 

J.  R.  Airey  and  I.  W.  M.  Smith,  J.  Chem.  Phys.  57,  1669  (1972). 


33 


5 


.  P.  R.  Poole  and  I.  W.  Smith,  J.  Chem.  Soc.  Faraday  II,  73,  1434  (1977); 
1447  (1977). 

6.  K.  G.  Anlauf,  P.  H.  Dawson,  and  J.  A.  Herman,  J.  Chem.  Phya.  58,  5354 
(1973). 

7.  M.  A.  Kwok  and  R.  L.  Wilkins,  J.  Chem.  Phys.  60,  2189  (1974). 

8.  M.  A.  Kwok  and  N.  Cohen,  J.  Chem.  Phys.  61,  5221  (1974). 

9.  J.  A.  Blauer  and  W.  C.  Solomon,  Int.  J.  Chem.  Klnet.  _5,  553  (1973). 

10.  R.  M.  Osgood,  Jr.,  P.  B.  Sackett  and  A.  Javan,  Appl.  Phys.  Lett.  20,  469 
(1972);  _22,  254  (1973);  also  J.  Chem.  Phys.  JSO,  1464  (1974). 

11.  J.  J.  Hlnchen  and  R.  H.  Hobbs,  J.  Chem.  Phys.  63,  353  (1975). 

12.  J.  F.  Bott,  J.  Chem.  Phys.  ^5_,  4239  (1976). 

13.  J.  F.  Bott  and  R.  F.  Heldner  III,  J.  Chem.  Phys.  68,  1708  (1978). 

14.  R.  L.  Wilkins  and  M.  A.  Kwok,  J.  Chem.  Phys.  70,  1705  (1979). 

15.  R.  L.  Wilkins,  J.  Chem.  Phys.  7£,  2700  (1979). 

16.  R.  L.  Wilkins,  "Vibrational  Relaxation  in  the  HF-DF  and  DF-HF  Systems," 
J.  Chem.  Phys.  (to  be  published). 

17.  R.  L.  Wilkins,  J.  Chem.  Phys.  5838  (1977). 

18.  E.  B.  Turner,  G.  Emanuel,  and  R.  L.  Wilkins,  The  NEST  Chemistry  Computer 
Program,  TR-0059(6240-40)-l ,  The  Aerospace  Corporation,  El  Segundo, 
California  (30  June  1970). 

19.  H.  K.  Shin  and  V.  H.  Kim,  J.  Chem.  Phys.  _64^  3634  (1976). 

20.  G.  D.  Billing  and  L.  L.  Poulsen,  J.  Chem.  Phys.  68,  5128  (1978). 


21.  C.  Clendening,  J.  L.  Stelnfeld,  and  L.  E.  Wilson,  Information  Theory 


Analysis  of  Deactivation  Rates  in  Chemical  Lasers,  AFWL-TR-76-144,  Air 


Force  Weapons  Laboratory,  Kirtland  AFB,  N.  Mex.  (October  1976). 


34 


22.  M.  A.  Kwok,  private  communication. 

23.  J.  F.  Bott,  Chem.  Phys.  Lett.  23,  335  (1973). 

24.  K.  Ernst,  R.  M.  Osgood,  A.  Javan,  and  P.  B.  Sackett,  Chem.  Phys.  Lett. 
23,  553  (1973). 

25.  J.  F.  Bott,  J.  Chem.  Phys.  70,  4123  (1979). 


/ 


35 


Appendix  A.  Rate  Coefficients  for  (R  >  R",  T"  ),  v  v'  ,  and  v  ♦  R 
Energy  Transfer  in  HF(v.  <  4)  +  HF(v2  *  0)  Collisions  at 
T  =  300  K  1 


Rotational  Relaxation  (R  -*■  R",  T') 


f 

HFCvj  Jj)  +  M  HFCvj,  Jj  -1)  +  M 


li 

il 

kf,  10^cm^/(mole  sec) 

0-4 

1-5 

20.0 

0-3 

6-10 

10.0 

0-3 

11-15 

5.0 

0-2 

16-20 

2.5 

0-1 

21-25 

1.25 

0 

26-30 

0.625 
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Vibrational  Relaxation  (v  ♦  R) 


HFCvp  Jt)  +  +  HF<v',  J'  )  +  Mx 


kf,  10l0cm"V(mole  sec) 


10-16 

10-16 

10-16 

10-16 

10-16 

10-16 

10-16 

10-16 

10-16 

10-16 

10-16 

10-16 

16-21 

16-21 

16-21 

16-21 

16-21 

16-21 

10-15 

10-15 

10-15 

10-15 

10-15 

10-15 

16-20 

16-20 

16-20 

16-20 

16-20 

16-20 

21-26 

21-26 

21-26 

21-26 

21-26 

21-26 

10-15 

10-15 


15.0 

38.0 

120.0 

52.0 

30.0 

30.0 


Vibrational  Relaxation  (v  ■*•  R)  (Continued) 


HF(Vl,  JL)  +  Ml  v 

V1 

J1 

vi 

4 

2 

3 

10-15 

4 

3 

3 

10-15 

4 

4 

3 

10-15 

4 

5 

3 

10-15 

4 

0 

2 

15-20 

4 

1 

2 

15-20 

4 

2 

2 

15-20 

4 

3 

2 

15-20 

4 

4 

2 

15-20 

4 

5 

2 

15-20 

4 

0 

1 

20-25 

4 

1 

1 

20-25 

4 

2 

1 

20-25 

4 

3 

1 

20-25 

4 

4 

1 

20-25 

4 

5 

1 

20-25 

4 

0 

0 

25-30 

4 

1 

0 

25-30 

4 

2 

0 

25-30 

4 

3 

0 

25-30 

4 

4 

0 

25-30 

A 

s 

n 

9s_,»n 

HF(v',  J')  +  Mj 

k^,  10^cm^/(mole  sec) 


70.0 

49.0 

28.0 

28.0 

5.8 

29.0 

58.0 

41.0 

23.0 

23.0 

3.7 

19.0 

37.0 

26.0 

14.0 

14.0 

4.1 
21.0 
41.0 
29.0 
16.0 
i  a  n 
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Appendix  B.  Rate  Coefficients  for  (R  -*•  R',  T'  ),  v  ♦  v' ,  and  v  R 

Energy  Transfer  in  DF(v^  <  3)  +  DF(v2  ”  0)  Collisions  at 
T  =  300  K 


Rotational  Relaxation  (R  +  R",  T'  ) 
kf 

DFCvj,  Jj)  +  M  DF  ( v  j ,  J:  -1)  +  M 


V1 

J1 

kf,  10^cm^/(mole  sec) 

>-4 

1-5 

15.0 

1-3 

6-10 

7.5 

1-3 

11-15 

3.75 

>-2 

16-20 

1.875 

1-1 

21-25 

0.9375 

0 

26-30 

0.46875 

Vibrational  Relaxation  (v  +  v') 

k, 

'1’  J  ] 

i)  +  DF(v2  =  0,  J,)  + 

DF(Vl  -1,  Jj)  +  DF(v2  =  1,  JL) 

k. 

d 

J1 

kf,  10  ^cm^/  (mole 

sec)  k^,  10^cm^/(mole  sec) 

0-5 

10.6 

13.2 

0-5 

6.5 

8.6 

0-5 

4.1 

5.5 
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